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SELF-CONSISTENT PHYSICAL PROPERTIES OF CARBON NANOTUBES 

IN COMPOSITE MATERIALS 


R. B. PEPES 1 , S J.V. FRANKLAND 2 ’ 3 , P. HUBERT 4 , AND E. SAETHER 5 

Abstract. A set of relationships is developed for selected physical properties of single-walled 
carbon nanotubes (S WCN) and their hexagonal arrays as a function of nanotube size in terms of the chiral 
vector integer pair, (n,m). Properties include density, principal Young’s modulus, and specific Young s 
modulus. Relationships between weight fraction and volume fraction of SWCN and their arrays are 
developed for polymeric mixtures. 

Key words, carbon nanotubes, density, modulus, hexagonal array, volume fraction, weight 
fraction, composites 

Subject classification. Materials Science 

Nomenclature 

ai unit vector 

a 2 unit vector 

b C-C bond length 

E a Young’s modulus of a SWCN (single-walled carbon nanotube) 

Ena Young’s modulus of a hexagonal array of SWCN 
E pa SWCN array specific modulus 

E pa SWCN specific modulus 

M w atomic weight of carbon 

m number of carbon atoms in a 2 direction 

n number of carbon atoms in ai direction 

N number of carbon atoms per unit length 
N a A vogadro ’ s number 

R n SWCN radius 

Rna effective SWCN radius in and hexagonal array 

R ne effective SWCN radius 

V a SWCN volume fraction of the hexagonal array 

V n SWCN volume fraction 

W a SWCN array weight fraction 

W n SWCN weight fraction 

Y graphene Y oung’ s modulus 
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A chiral vector function 

A. equilibrium separation distance between SWCN 

v equilibrium separation distance between S W CN and polymer 

p a SWCN array density 

p m SWCN/polymer composite density 

p n SWCN density 

Poa SWCN density in a hexagonal array 

p p polymer density 

p a SWCN effective array specific gravity 

p n SWCN effective specific gravity 

1. Introduction. The use of single-walled carbon nanotubes (SWCN) to further the reinforcement 
of polymeric materials is one of the current applications being explored [1,2] for this new material form 
since its discovery a decade ago [3]. In determining the nanotube contribution to the overall properties of 
a polymer composite, it is important to note that the ‘effective fiber’ is the nanotube or nanotube array, 
including their enclosed volume. This basic consideration affects a number physical properties of the 
composite such as the density, modulus, specific modulus, and the conversion between weight and 
volume fraction. 

Considerable attention has been given to measuring and modeling the principal Young’s modulus 
of carbon nanotubes. Some of these studies assume a nanotube shell, and others the enclosed nanotube 
volume. These different geometry definitions have led to reported values differing by almost two orders 
of magnitude [4-9]. For example, an experiment, determining the Young’s modulus of the SWCN from 
vibrational analysis, assumed in the data analysis a nanotube shell of thickness 0.34 nm resulting in a 
Young’s modulus of 1.25 TPa for SWCN of 1.0- 1.5 nm in diameter [6]. . AFM (atomic force 

microscopy) deflection experiments, which included the entire cross-section of the nanotube in the 
analysis, measured a Young’s modulus of 1.31 TPa for an array 3 nm in diameter with the modulus 
dropping to 67 GPa for an array 20 nm in diameter [7]. The average diameter of individual SWCN in this 
study was 1.4 nm. The geometrical dependence of the nanotube modulus has been confirmed with 
theoretical calculations. Non-orthogonal tight-binding calculations based on a nanotube shell at a 
thickness of 0.34 nm, but parametrized for any shell thickness, yield 1.22-1.26 TPa for SWCN of 0. 8-2.0 
nmin diameter [8]. When the entire nanotube cross-section is included, molecular dynamics simulations 
of nanotubes in compression that utilize the Tersoff-Brenner potential for the C-C bonds, yield a Young s 
modulus of 1.5 TPa at SWCN diameter of -0.5 nm which reduces to 0.2 TPa at a SWCN diameter of 4.0 
nm [9]. 

The density of carbon nanotubes has been measured experimentally for randomly oriented 
nanotube arrays under high pressure [10]. Extrapolating to 1 atm, the density is close that expected of a 
perfect crystal lattice of SWCN reported as 1.33 g/cm 3 . Another study reports a calculation of nanotube 
density based on assuming a shell geometry [11]. In the present work, a simple expression is derived for 
the density of an individual SWCN and a perfectly hexagonal crystalline array of SWCN in terms of both 
the radius and the chiral vector. With a given expression for density, it is further possible to develop a 
conversion between weight fraction and volume fraction of SWCN in a composite. This conversion is 
especially useful as a connection between theoretical predictions which tend to be in volume fraction, and 
experimental measurements which tend to be in weight (or mass) fraction. 

One last consideration in the properties that will be reported here is the definition of nanotube 
size in the composite. Due to the chemical interactions of the polymer and the nanotubes, there is a small 


equilibrium separation distance between the nanotube and either the surrounding polymer or other 
nanotubes in the SWCN array. This equilibrium separation distance is especially significant because the 
diameter of the nanotube is on the order of nanometers across. In this paper, the equilibrium separation 
distance between the nano tube and polymer is arbitrarily taken as O.o4 nm, the graphitic plane separation 
distance [12,13]. A value calculated in this work is used for the nanotube-nanotube equilibrium 
separation. 

In this work, analytical expressions are derived for nanotube density, modulus, specific modulus 
and weight-volume fraction conversion assuming that the geometry of the nanotube or nanotube arrays 
includes all of the enclosed volume. This volume is dependent on the size of nanotube which is in turn 
dependent on the chirality (helicity) of the nanotube, generally described by the indices of its chiral vector 
(n,m) [14]. After a general description of nanotube geometry and nomenclature, each property will be 
addressed in turn for both SWCN and SWCN arrays: density, modulus, specific modulus, and weight- 
volume fraction conversion. These formulae provide a consistent set of relationships for use by 
researchers interested in calculating the effects of nanotube size on the effective properties of SWCN 
composites. 

2. Geometry. 

2.1. SWCN. The SWCN has been described as a single graphene sheet rolled up with varying 
degrees of twist as described by its chiral vector, C h : 

(2.1) C h =na x fma 2 

where ai and a 2 are the unit vectors of the two-dimensional hexagonal lattice and the indices ‘V’ and "m” 
are integers [14]. These vectors and their summation axe illustrated in Figure 2.1 for 3 specific examples 
of nanotubes. The solid lines perpendicular to the vector sum of a : and a 2 form the seam of the wrapped 
graphene sheet. The chiral vector C h is also referred to by its indices as (n,m). 

Nanotubes with chiral vectors of (n,n) and (n,0) have no twist and are classified as achiral 
nanotubes. These two special cases are sometimes denoted armchair and zig-zag respectively, referring 
to the pattern of the carbon atoms around the nanotube circumference [14]. Specific examples are 
included in Figure 2.1(a) for the vector diagram and structure (front and side views) of (6,6) armchair 
nanotube, and in Figure 2.1(b), for the (10,0) zig-zag nanotube where n=10. SWCN that do not fit either 
of these classifications are referred to as chiral nanotubes [14]. A (6,5) chiral nanotube is shown in 
Figure 2.1(c). 

The SWCN radius, R n , is a function of the integer pair, n and m, and the C-C bond length b [14]: 

(2.2) R n =^-j3 A 


where 



For the SWCN, the C-C bond length is equal to 0. 142 nm [14]. This radius, does not the equilibrium 
separation distance, v, between the nanotube and the surrounding polymer. Therefore, an effective radius, 
R ne is defined to include v as 

(2.4) R ne =^-y[3A + ^~ . 

l7t 2 

The relation of R n> R ne> and v y is shown in Figure 2.2(a) for the SWCN, and the resulting effective fiber 
cross-section in Figure 2.2(b). The surrounding polymer and the nanotube are each assigned one half of 
the equilibrium separation distance. The value of v is dependent on the properties of the particular 
polymer, especially which functional groups of the polymer repeat unit are in the vicinity of the nanotube 
and how they align with respect to the nanotube. In the calculations a default value of 0.34 nm is used 
which is the separation of graphite planes [12, 13]. Other values from simulation force fields are 
considered. 



a) Open cylinder model of SWCN b) Filled cylinder model of SWCN 


FIGURE 2.2. SWCN” nomenclature. 


2.2. SWCN Hexagonal Array. The synthesis of SWCN typically generates arrays of SWCN in 
a hexagonal cross-sectional arrangement [15]. For a SWCN in an array, the effective radius is 

b X 

(2.5) R na =—^A + — 

2n 2 

where X is the nanotube-nanotube equilibrium separation distance within the array analogous to V 
specifically defined, in this work, for the polymer. The relation of R n , Rna, and >4 is shown in Figure 2.3(a) 
for the SWCN array, and the resulting effective fiber cross-section in Figure 2.3(b). 

The value of X was calculated for various SWCN arrays from molecular statics simulations using 
the method in Ref. 16. In these simulations, the potential energy of the array was minimized, using the 
Lennard-Jones potential (with parameters of 8=34.0 K and a=0.3406 nm) to describe the nanotube- 
nanotube interactions. The results, summarized in Table 2.1, suggest that the separation distance is not 
strongly dependent on SWCN diameter, but rather can be taken as a constant equal to 0.318 nm in 
agreement with the equilibrium distance of 0.313 nm determined in Ref. [17] and 0.315 nm in Ref. [22] 




Table 2.1 

SWCN (n,0) Hexagonal Array Separation Distance. 


n 

Diameter 

(nm) 

Separation 
Distance, X 
(nm) 

6 

0.48 

0.316 

12 

0.94 

0.317 

18 

1.41 

0.317 

24 

1.88 

0.318 

54 

4.23 

0.318 

96 

7.51 

0.319 


for a similar range of SWCN diameters. For a constant equilibrium separation distance, the SWCN 
volume fraction of the hexagonal array, V lt is a constant equal to the volume packing fraction of 0.906: 

(2.6) V a =-^j= = 0.906 




a) Carbon nanotube array 


b) Effective reinforcement array 


FIGURE 2.3. SWCN array nomenclature. 
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3. Density. The density of SWCN is here derived from the mass of carbon atoms in the 
nanotube lattice per unit volume enclosed by the nanotube cylinder shown in Figure 2.2(b). The number 
of carbons per unit length is 


(3.1) 


N = 


4A 
3 b 


Combining equations (2.4) and (3.1), the nanotube density, p Di can be expressed in terms of the chiral 
vector as follows: 


NM W _ \6tzAM w 

^ R l N a ~ 3b{ib 2 A 2 + l-Bbntv v + ttV )N a 


where M w is the molecular weight of carbon, N a is Avogadro’s number, and R ne the effective radius which 
includes v. For any particular polymer composite, the SWCN density is therefore dependent on the radius 
of the nanotube. This dependence of SWCN density on diameter is plotted in Figure 3.1 for v=0.34 nm. 
Densities of specific nanotubes are marked with solid symbols. Over a range in diameter of 1-14 nm, 
that SWCN density decreases with increase in diameter by an order of magnitude. The SWCN density 
contribution is also dependent on the separation distance v. In Table 3.1, a specific example of the 
sensitivity to v is given for a (10,10) SWCN. For reference, the equilibrium separation distance for 
crystalline and amorphous polyethylene-nanotube composites with and without chemical bonds between 
the polymer and the nanotube is 0.38-0.42 nm. This range comes from molecular dynamics simulations 
of these composites and is determined from the separation of the nanotube peak from the first polymer 
peak in the radial distribution function. Details of the simulations can be found in Ref. [18]. The starting 
value of 0.25 nm was chosen because it is a typical hydrogen-hydrogen sigma parameter for the Lennard- 
Jones potential [19]. 


Table 3.1. 

Density of a (10,10) SWCN in Different Matrices . 


v nm 

Density in g/cm 3 

0.25 

1.602 

0.30 

1.506 

0.32 

1.471 

0.33 

1.453 

0.34 

1.436 

0.35 

1.419 

0.36 

1.403 
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Figure 3. 1 . SWCN density versus diameter. 

The density of a SWCN array is calculated by dividing the density of a single SWCN by the 
volume fraction of SWCN in the array which includes effective volume increase due to the nanotube- 
nanotube equilibrium separation distance X. From equations (2.5), (2,6), and (3.1), SWCN array density 
is 


\6nh M w 

3&(3£ 2 A 2 + 7^3b7tAA + ]n o _ 

Like the density of an individual SWCN, the array density shows a significant reduction with SWCN 
diameter. The density of a SWCN array differs from that of a SWCN only by a factor of the maximum 
volume fraction (0.906) and the difference in nanotube-polymer and nanotube-nanotube equilibrium 
separation distances. 

4. Principal Young’s Modulus. The principal Young’s modulus is taken in the direction of the 
longitudinal axis of SWCN and its arrays. It is assumed in the following that the nanotubes are 
continuous and that the arrays consist of SWCN of identical diameter. 

4.1. SWCN. A simple approach to calculating the Young’s modulus of a SWCN is to map the 
stiffness of a graphene sheet onto the enclosing surface of the nanotube cylinder. The concept of using 
the graphene stiffness as a key parameter in calculating the Young’s modulus of SWCN is not unique and 
has been used by others. (See for example, Refs. 7 and 20.) 


(3.3) 


Pa =VaP„a = V a 


NM 


n 


nRlN a 2V3 
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To derive an expression which relates Young’s modulus of the nanotube in terms of its chiral 
vector (n,m) as A directly to the graphene data, the modulus-area product of the effective modulus of the 
SWCN in the composite, E n , as a filled cylinder of radius R ne is set equal to that of a thin-walled cylinder 
of outer radius R ne , thickness v, and modulus Y of the graphene sheet. The effective modulus E n can then 
be written as 

MRS 

^ n 7 7 

4 +4R n v + v 


(4.1) 


where Y = 


(cfi-c | 2 2 ) 


c 
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where Cn and C12 are the stiffness constants of graphene [21]. Using equation (2.2) for R n , equation (4.1) 
is rearranged as 


(4.2) 


443tayav 

3b 2 A 2 +2j3b7rAv + x 2 v 2 


Equation (4.2) can be derived from equation (1) in Ref. 20 if Y is used rather than Cn. In Figure 4.1, the 
Young’s modulus of the SWCN from equation (4.2) is plotted as a function of diameter, taking Y=1029 
GPa as the value of the Young’s modulus of the graphene sheet and the separation distance as 0.342 nm 
[12]. Specific nanotube chiralities are marked with solid symbols in Figure 4.1. Equation (4.1) shows 
that the modulus of the SWCN decreases with increase in radius. A similar reduction in modulus with 
nanotube radius was derived from molecular dynamics simulation results of SWCN moduli using the 
Tersoff-Brenner potential to describe the C-C bonds [9]. In Table 4.1 a comparison of specific nanotube 
chiralities from equation (4.2) with the results of the graphene sheet model in Ref. 20 and the expression 
derived from molecular dynamics simulation results in Ref. 9 shows that the all three models provide very 
similar results. 

Table 4.1 

SWCN Young ’s Modulus for Specific Nanotube Chiralities. 


Chiral 

n 

Indices 

m 

1 

E„ 

Eq. (4.2) 
(GPa) 

SWCN 

E„ 

Ref. 9 
(GPa) 

E n 

Ref.20 

(GPa) 

Arrays 

E a 

Eq.(4.3) 

(GPa) 

E a 

Ref. 7 
(GPa) 

10 

10 

662 

642 

680 

574 

608 

18 

0 

647 

618 

664 

561 

594 

24 

0 

536 

465 

550 

462 

492 

50 

0 

304 

228 

311 

259 

279 

96 

0 

171 

123 

175 

145 

157 
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FIGURE 4.1. Dependence ofSWCN Young ’s Modulus on nanotube size . 

4.2 SWCN Hexagonal Array. The Young’s modulus, E a , of a SWCN array containing SWCN 
with moduli described by equation (4.2) is obtained from equations (2.5), (2.6), and (4.2) as 

2n 2 bYM 

(43) “ ' “ “ " 3i> 2 A 2 + iSbTTKX + 7PX 1 


Values for the principal modulus of the SWCN array shown in Figure 4.2 are calculated from equation 
(4.3) and show that the array moduli follow identical trends to the individual SWCN. The array moduli 
differ from the SWCN moduli only by the product of the maximum volume fraction (0.906) and the 
difference in van der Waals distances (A=0.318 nm vs. v = 0.34 nm). In Table 4.1, a comparison is 
included of moduli calculated from equation (4.3) with those from the array graphene sheet model for 
SWCN arrays in Reference [7]. The model in Reference [7] differs from equation (4.3) by using C n 
instead of Y and by using separate values for the graphene sheet thickness and the nanotube-nanotube 
separation distance. Furthermore, the values from equation (4.3) show excellent agreement with those 
calculated via a lattice dynamics method for selected armchair (3n,3n) and zig-zag (3n,0) SWCN [22]. 
Selected data points (solid symbols) from Ref. 22 are plotted with equation (4.3) in Figure 4.2. 

5. Specific Axial Modulus. The specific modulus, Zt prli of the SWCN is defined as the 
ratio of the principal Young’s modulus, to the dimensionless magnitude of the density, p n , analogous 

to specific gravity. Combining equations (3.2) and (4.2) yields 
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(5.1) 


E 


pn 


E n _ 3'j3b 2 YvN a 
K~ 4M w 


Since the same geometry is utilized for calculating both modulus and density, the resulting equation (5.1) 
yields results independent of the SWCN diameter. 

In a similar manner the specific modulus for the SWCN array, E ^ , is obtained by combining 
equations (3.3) and (4.4): 


(5.2) 


Pa 


3y5 b 2 YAN a 
P a 4 M w 



FIGURE 4 . 2 . Dependence of SWCN array axial modulus on nanotube size. 


These results show that the specific modulus of the SWCN array is, with the exception of the separation 
distance parameter, identical in form to that of the SWCN. 

6. Weight Fraction/Volume Fraction Conversion. With an expression for the density of the 
SWCN, a relationship can be derived between weight fraction and volume fraction of SWCN in a 
composite. Consider a SWCN/polymer mixture of density, p m with a SWCN volume fraction, V a , a 
SWCN density of p n and a polymer density of p r The SWCN volume fraction can be expressed as 
follows: 
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( 6 . 1 ) 


Pm-Pp 

P n ~Pp 


It is also possible to express the SWCN volume fraction in terms of its weight fraction, W m by definition: 


( 6 . 2 ) 


w a 

P„ 


w„p p 

W n p p +{1-W n )p n 


In terms of the chiral vector, equation (6.2) becomes 


(6.3) 


W.3M/J 

3b 2 A 2 +2^[3bzAv + 7c 2 v 2 ) 

'Pp 

W nPp 3bN a \ 

3b 2 A 2 

\ 

+ 2y[3b7rAv + 7T 2 V 2 j 

1+ (l—W n )16 tvAM w 


This relationship between volume fraction and weight fraction for SWCN-polymer mixtures is plotted in 
Figure 6.1 (solid lines) for (6,6), (12,12) and (18,18) SWCN for a polymer density of 1 g/cm 3 in 
comparison to the V n =W n line (dotted). These curves show that as the diameter of the SWCN is increased 
the non-linearity in the relationship decreases. The influence of variations in polymer density upon the 
relationship between volume and weight fraction is shown in Figure 6.2 for the (10,10) SWCN in 
polymers of 0.8, 1.0, and 1.2 g/cm 3 . Here increases in polymer density also reduce the non-linearity of 
the curves. 

For SWCN arrays in a polymer mixture, a similar relationship between weight-fraction, W a , and 
volume-fraction V a can be derived as 

W a 3^3 bN a (3b 2 A 2 + 2j3bn . A v + n 2 X 2 )p p 
^ 6 ' 4) V “ ~ W nPp 3j3bN a [3b 2 A 2 + 2j3bizA v + n 2 X 2 )+{l -W n )8x 2 AM w 

8, Conclusions. The overall set of analytical expressions presented here for the density, moduli, 
and weight-volume fraction conversion of nanotubes and nanotube arrays are derived as functions of the 
nanotube radius or chiral vector (n,m), and includes the equilbrium separation distance in the nanotube 
geometry to obtain the effective enclosed volume of the individual nanotube or array. The density and 
moduli of SWCN in polymer and SWCN arrays differ only by the array packing fraction and the 
equilibrium separation distances due to the chemical interactions between the nanotube and polymer or 
other nanotubes. 
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, Volume Fraction V n , Volume Fraction 



Figure 6.1. SWCN volume fraction versus weight fraction (polymer density = 1 g/cm 3 ). 



0 0.2 0.4 0.6 0.8 1 

W„ , Weight Fraction 


FIGURE 6 . 2 . SWCN (10,10) results for different polymer densities. 
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A simple model for the SWCN principal modulus is derived assuming that the nanotube consists 
of a hollow cylinder with properties of the cylinder wall equal to the stiffness of the graphene sheet. This 
model agrees well with other graphene sheet models for SWCN and arrays in the literature, as well as 
published molecular dynamics simulation and lattice dynamics results. The modulus calculated from this 
model for both the SWCN and array shows an order of magnitude decrease over an increase in diameter 
from 1 to 14 nm. With this model, the specific modulus of both the SWCN and its hexagonal arrays 
was found to be independent of SWCN diameter. 

Finally, a set of equations was derived for the relationship between weight fraction and volume 
fraction for SWCN and hexagonal arrays. The relationship requires knowledge only of the components of 
the chiral vector or radius of the SWCN and the weight fractions of the constituents to determine volume 
fraction of the SWCN or its hexagonal arrays in a mixture. For larger nanotubes in the same medium or 
the same size nanotubes in denser media the conversion curves becomes more linear and closer to a 1:1 
conversion, assuming idealized, undeformed nanotubes. 
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